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Abstract 

In a recent study of large non-null sample covariance matrices, a new sequence of functions gener- 
alizing the GUE Tracy- Widom distribution of random matrix theory was obtained. This paper derives 
Painleve formulas of these functions and use them to prove that they are indeed distribution functions. 
Applications of these new distribution functions to last passage percolation, queues in tandem and to- 
tally asymmetric simple exclusion process are also discussed. As a part of the proof, a representation of 
orthogonal polynomials on the unit circle in terms of an operator on a discrete set is presented. 



1 Introduction 

Let Ai(u) denote the Airy function. It has an integral representation 

Ai(u) = — I e^ ua+ ^da (1.1) 
2n J 

where the integral is over a curve from ooe 5i7r / 6 to ooe i7r / 6 . The Airy kernel (see, e.g. ^JEZj) is defined as 
A/ . M(u)M'(v)-M(u)M'(v) f°° ,., sw 

A(u,v):=— — — — — =/ Ai(u + z) M(z + v)dz. (1.2) 

u — v J 

Let A x be the Airy operator acting on L 2 ((x, oo)) whose kernel is given by A(u, v). Define 

F (x) :=det(l-A K ). (1.3) 

The 'GUE Tracy- Widom distribution function' Fq(x) is the limiting distribution function of various models 
in mathematical physics, probability and statistics (see e.g. and references in it). 1 Especially in statis- 
tics, the largest eigenvalue of the sample covariance matrix of complex Gaussian samples with the identity 
covariance (the so-called null case) is known to have the limiting distribution given by Fq(x). An intriguing 
result by Tracy and Widom |27j is that the Fredholm determinant has an alternative expression: 



F {x) = det(l - A*) = exp(- J 



(s-x)u z (s)ds), (1.4) 
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1 In many literatures, Fq is denoted by F2. In this paper, we reserve F2 for a different function. 



where u{x) is the solution to the Painleve II equation 

u" = 2u 3 + xu, (1.5) 

subject to the condition 

u(x) ~ — Ai(x) as x — > +00. (1.6) 

It is known that there is a unique global solution to the equation (|1.5|) with the condition (|1.6|) , and the 
solution satisfies (see, e.g. jl 71 IIP) 1 ) 

/e-^ 3/2 \ 

u(a;) = -Ai(x)+0\ /4 j, a; -> +00 (1.7) 

«(*) = ~\/t^( 1 + (^))' 21 ^ _o °- ( L8 ) 

Recall that Ai(x) ~ e~* x 1 /(2v / 7rx 1 / 4 ) as x — > +00. The right-hand-side of (|1.4(l provides a practical 
formula to plot the graph of Fq numerically. 

For m = 1, 2, 3, . . . and for complex numbers w\, 11)2, ■ ■ ■ , define 

1 f m 

s^(u;wi,...,w m )=s^ m \w u ...,w m ):=— e ¥- 3 +^T\ da ( L9 ) 

Z7T / -*• x k;, + ia 

where the contour is from ooe 5 **"/ 6 to ooe 171 "/ 6 such that the poles a — iwi, . . . ,iw m lie above the contour. 
Also define 

, « m— 1 

t^ m \v;w 1 ,...,w m - 1 )=^ m \w 1 ,...,w m . 1 ):=— / e ¥" 3 +^ JJ (w , _ ib)db (1.10) 

^ ^ i=i 

where the contour is from ooe 5l7r / 6 to ooe I7r / 6 . Comparing with Ijl.ip . v- m ^ is a sum of derivatives of the 
Airy function. On the other hand, when w\ = ■ ■ ■ = w m = 0, s'" 1 * 1 is a sum of anti-derivatives of the Airy 
function. However for general lUj's, s^" 1 ) is a Cauchy-type transform of the integrand of the Airy function. 
Define 

Fk(x;wi, . . . , Wk) 

( 1 ( 1 n \ (1-11) 

:= Fq(x) ■ dct \5mn- < _ t s w (wi, . . . , ttf m ),t w (u;i, . . . ,io„_i) >t2(( I;OC )) I 

V 1 A 2 / i< m ,n<k 

where <,>(j; ; oo) denotes the real inner product in L 2 ((x, 00)); 

< S (m) (wi, ■ • . , W m ),t (n) (wi, . . . (ttfn-l) >Z,2((x,oo)) 

! x (1-12) 

s (m) (wi, . . . , w m ) I (ti)i (n) (it; wi, ■ ■ .,w n -i)du. 



1 - A : 

(It is well-known that 1 — A x is invcrtible.) Set 

F k {x):=F k {x;0,0...,0), fc = l,2,.... (1.13) 
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The functions Fk(x; Wx, . . . , Wk) were introduced recently in £Q as limits of the distribution functions of the 
largest eigenvalues of certain non-null complex sample covariance matrices and also other probability models. 
Sec Section [21 below for more details on the motivations. The purpose of this paper is to find a Painleve 
type formula for Fk(x; w\, ■ ■ • , Wk) analogous to l|1.4|) . Such formula is used to prove that Fk(x; w\, ■ ■ ■ , Wk) 
is indeed a distribution. It also allows us to be able to plot the graph of Fk(x). 

1.1 Results 

1.1.1 Alternative determinantal formula 

We first obtain an alternative determinantal formula of Fk(x; Wx, . . . , lUfc). The definition involves the 

functions and t^" 1 ' and it is not transparent that the formula is symmetric in w\, ■ . ■ , wj~, which should 
be the case from its origin in the sample covariance matrix pQ (see also Section [2 below). This symmetry is 
clear in the following theorem. 
For a complex number w, set 

C w (u) := i / e^+^—^—da (1.14) 
Z7T J w + ia 

where the contour is, as in the definition l|1.9f) of s^ m \ from ooe 5l7r / 6 to ooe l7r / 6 such that the pole a = iw 
lies above the contour. Hence s^(u;wi) = C Wl (u). Also note that t^\v) = Ai(v). 

Theorem 1.1. With above notations, for real x and complex w set 

f(x,w) := 1- < - ^ C w ,Ai >L2((x,oo))= 1- J ( 1 ^ A Cw\ (u) Ai(u)du. (1.15) 
For distinct complex numbers wx, . . . , Wk, 



det 



[(w m + D x ) n - 1 f{x,w m )y 



F k {x;wx,...,w k ) =Fo{x) — - ■ = (i.ie) 

(w n - w m ) 

l<m<n</c 

where D x = ^ denotes the derivative with respect to x. When some of Wj 's coincide, the above formula still 
holds by using the I'Hosptial's rule for the right-hand- side of (|3.1() . 

Remark. P. Deift and A. Its pointed out that this formula resembles the Darboux transformation in the 
theory of integrablc systems (see e.g., |23|). It would be interesting to identify the above formula in terms 
of a Darboux transformation of an integrablc system. 

This theorem follows from row and column operations of 1)1.11(1 exploiting the fact that t^ n > is a sum of 
derivatives of Ai and that s( m > is a linear combination of C Wj . The proof is given in Sectional 

1.1.2 Painleve formula 

In the next theorem, we show that the function f(x, w) defined in (|1.15fl is related to the Painleve II equation. 
First we need a definition. Let M(z; x) = ( ^!(z) M^zj ) ^ e ^ ne 2x2 matrix-valued solution to the 



following Riemann-Hilbert problem: 
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• M(z; x) is analytic for z € C \ R and is continuous for z 6 C \ ' 

• For z G R, 

M+fe x) = M_(z;x) f s s ~ 6 (1.17) 

where M + (z;x) (resp. M_(z;x)) denotes the limit of M(z';x) as z' — » z from the bottom (resp. top) 
of the contour R. 

• M(z; x) — > / as z — > oo. 

The precise statement of the last condition is the following: there is e > such that M(z;x) = I + 0(z~ 1 ) 
uniformly as z — > oo for z in sectors e < Arg(z) < n—e and 7r + e < Arg(z) < 2n— e, and M(z; x) is bounded 
for all z G C \ R. 

This is the Riemann-Hilbert problem for the Painleve II equation when the so-called monodromy data 
satisfies p = — q = 1 and r = [lHlllllfTHj . It is known that there is a unique solution to this Riemann-Hilbert 
problem. Moreover, as z — * oo, there is an expansion of form 

M { z;x) = I + ^ + o(% M lW = iH Ui f\ (1.18) 
z \z J 2i \-u{x) v(x)J 

where u(x) is the solution of the Painleve II equation (|1.5|) satisfying i|1.6fl . and 

v(x) = [ u{sfds. (1.19) 

J oo 

The following theorem shows that f(x,w) is expressible in terms of the Riemann-Hilbert problem for 
Painleve II equation. 

Theorem 1.2. The junction f(x,w) defined in (I1.15|) satisfies the following: 

(M 22 (-^w;x), Re(w)>0 
j{x,w) = < (1-20) 
[-M 2 i(-^iw;x)ei w3 - xw , Re{w) < 0. 

Note that from the jump condition i|1.17[l . f{x, w) is continuous for w G R, and hence is an entire function 
in w. 

Together with Theorem 11.11 Theorem 11.21 yields the desired Painleve II formula of Fk(x; wi, . . . , Wk), 
which is the main result of this paper. 

Corollary 1.3. The function Fk(x; w\, . . . , Wk) defined by (|1.11|) is equal to (|1.16fl with f{x,w) given 
by (JTUJ. 

The function given in the right-hand-side of (|1.2U|) had previously appeared in ^ (equation (2.22)) 
and [3j (equation (3.5)) as a limiting function for a last passage site percolation model. In the context of 
symmetrized random permutations and last passage percolation models, 01 0j showed, among other things, 
the k = 1 case of Corollary 1 1.31 

F 1 (x,w 1 )=F (x)f(x,w 1 ) (1.21) 

where f(x,w) given by the right-hand-side of i|1.20fl . This paper proves that the general case is expressible 
in terms of derivatives of the same function f(x, w). 
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1.1.3 Properties of f(x.w) 

The papers 0] |3] proved several properties of the function defined by the right-hand-side of (|1.2U|) . By 
setting w i— ► iw in Lemma 2.1 of 0] or Lemma 3.1 of we find the following properties of f(x,w). The 
following complementary function is useful: set 

{M 12 {-\iw;x), Re(w)>0 
g(x,w):=i (1.22) 

[-M n (-±iw;x)e3 w ~ xw , Re(w) < 0. 

Lemma 1.4 ( |3"| ) . The following holds, 
(i). f(x,w),g(x 7 w) are real for w G R. 
fnj. For eac/i /ia;ed w; G C, hs i' ^ +oo 



r 3/2 



and dsn — oo, 



/(*,«;) = l + 0(e~ c * ), (1.23) 
g{x,w) = -e^ w3 - xw (l + 0(e- cxa/2 )) (1.24) 



f(x,w) ~ _l e ^'MM 3/2 +lM™-» 2 M 1/2 (1.25) 
3 (x, W ) ~ __ L e *« 8 -*M 3/a +*W«-» 2 I-I 1/2 . (1.26) 



(Hi). 



V2 



lim /(x,ii;)=l, lim g(o;,u;)=0, (1-27) 

'UJ— > + oo W — >-+oc 



lim f(x,w)=0, lim g(x,w)=0, (1.28) 

10 — * — oo w — ► — oo 

f(x,0) = £(x), g(x,0) = -S(x), (1.29) 
E(x):=exp^J u(s)ds\. (1.30) 

(£(x) was denoted by E(x) 2 in 



(iv). For all x G K and w £ C, 











-1 




w)J 


- ( 






w)\ 


■ ( 


-1 




w)J 





u(x)\ H(x,w) 

u(x) -w J \g{x,w))J' 

{u{x)) 2 —wu(x) — u'(x)\ / f(x,w) 

-wu(x)+u'(x) w 2 - x - (u(x)) 2 ) \g(x,w) 



(v). 



(1.31) 
(1.32) 



f(x,w) = - 5 (x,-«;)e^ — , (1.33) 
fffoto) = -f(x,-w)ei w3 - xw . (1.34) 
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(vi). For each fixed as w — > — oo, 

i r v 



f(y^\w\+w 2 ,w) -> ed(y) = -y=l e~*' ds, 



OO 



g(yV\M+^ 2 ,w) ~ _ e iM 3 +v%M 3 / 2 . (1.36) 

Note that (|1.31|) and (|1.32|) are the Lax pair equations for the Painleve II equation. Hence Theorem ll.il 
Theorem II . 21 and Lemma 1 1 . 41 yields that Fk(x;w\, . . . ,Wk) is expressible in terms of the Lax pair equations 
of the Painleve II equation. 

Remark. After this paper was completed, Harold Widom found a different proof of l|1.31fl and (|1.32|) for 
f(x,w) defined by l|1.15f) and g(x,w) defined by g(x,w) = — ( jits C w )(x) using the method of j23|. The 
proof of Widom is algebraic and is more direct. On the other hand, the current paper proves a general 
identity and then takes a limit as outlined in subsection II . 21 below. 

From l|1.31[) and (|1.32f> . f(x, w) itself satisfies a second order linear differential equations in x and w with 
coefficients involving u{x). 

Corollary 1.5. Denoting J^f(x,w) = f'(x,w) and -§^f{x,w) = f(x,w), f satisfies 

-f"+(-~w)f' + u 2 f = (1.37) 



and 



-f+ —+w'-xf+[ — +« 4 + xu z - (u'Y)f = 0. (1.38) 



Together with the initial conditions f(x, 0) = £{x) and f{x,0) = (u 2 (x) + u'{x))£ (x), l|1.38|l may provide 
a numerical way to compute the function f(x, w), and hence F^{x; W\, . . . , Wk)- 

1.1.4 Formula of Fk(x) 

When w\ = ■ ■ ■ = Wk = 0, using the l'Hospitals' rule in (|1.16|) . 



F k (x) = F k (x; 0, 0, . . . , 0) = -1—F (x) ■ dct (d^^w + D^ffa w)}\ w J) . (1.39) 

Hj = oj' ^ / X<m,n<k 



By using (|1.31|) . (|1.32|) and (|1.29fl . one can in principle compute the determinant. The first three functions 
are 

F 1 {x) =F (x)£(x), 

F 2 (x) = F (x)£(x) 2 {l + u(x + 2u 2 + 2u')}, (1.40) 
F 3 (x) = F {x)£(x) 3 {l + 2u(x + 2u 2 + 2u') + i(u 2 - u'){x + 2u 2 + 2u) 2 }. 

Using the numerical evaluation of the Painleve solution u(x) which is available at the website of M. Prahofcr 
(http : //www-m5 . ma . turn . de/KPZ) , these formulas provide a convenient way to plot the graphs of F k . Figure^ 
is the graphs of the density function j^Fj(x) for j = 0, 1, 2, 3. Note that the function moves to the right as 
the index k increases. The numerical means and the standard deviations of Fh(x) are the following: 
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Figure 1: Graph of ^Fj(x), j = 0, 1, 2, 3 from the left to the right 





mean 


standard deviation 


F 


-1.771... 


0.90... 


*i 


-0.494. . . 


1.11... 


F 2 


0.543.. . 


1.18... 


F 3 


1.445.. . 


1.21... 



1.1.5 Fk are distribution functions 

The Painlev'e formula obtained above can be used to prove that Fk (x; Wi , . . . , Wk) is indeed a distribution 
function. 

Corollary 1.6. The function Fk(x; w\, . . . , Wk) is a distribution for real Wi, . . . , Wk- 

Proof. In p^, the function Fk(x;wi, . . . ,Wk) are shown to be continuous, non-decreasing and converges to 1 
asi-> +oo (see the paragraph after (25)). We need to show that Fk(x; Wi, . . . , Wk) — + as x — > — oo. From 
(|1.25|l . all entries of the determinant in both (|1.16(l and (|1.39|l are in absolute value less than or equal to 
Q e -c\x\ 1 f or some constants C, c > as x — > — oo. Also Fq(x) < Ce~ c \ x \ for some other constants C, c > 
(see e.g. (2.13) of 0]). Hence Fk(x;wi, . . . , Wk) = 0(e~ c \ x \ ) for some constant c > as x — ► — oo. □ 

1.2 Outline of the proof and orthogonal polynomials on the unit circle 

Theorem 11.11 is obtained by applying a sequence of row and column operations to the determinant H.llfl . 
This part is the main bulk of the paper and the proof is given in Section [3j 

The proof of Theorem 11.21 is indirect. We use a representation of orthogonal polynomials on the unit 
circle in terms of an operator on a discrete set. Since such a representation may be interesting in itself, we 
state it here. This formula follows from a general identity (see H4.ll below) between Toeplitz determinants 
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and Frcdholm determinants on integer lattice obtained by Geronimo and Case |14| . and also independently 
by Borodin and Okounkov [7] (see also jHJ|S] for shorter proofs). 

Set No := N U {0}. Let cf>(z) be a function which is positive on the unit circle. For simplicity of 
argument, we assume that 4>{z) is analytic in neighborhood of the unit circle. Let <f> — 4 > +4>- be a Wiener- 
Hopf factorization of 4> where 4>+ extends to a non-vanishing analytic function interior of the circle and <f>- 
extends to a non-vanishing analytic function exterior of the circle. Set 

(, 41) 

For a function f(z) on the unit circle, fk denotes its k th Fourier coefficient: 

fk ■= / z- k f(z)--. (1.42) 

J\z\ = l 2mz 

Proposition 1.7. Let 7r n (z) be the monic orthogonal polynomial on the unit circle with respect to the measure 
0(z) Ijiriz ; an d ^ ^ni 2 ) = zn ^n(k) ^ e ^ s ^-transform. For <f> satisfying above conditions, 

<{z) = e-ZZLA^^L < _ _ _J__ Pn Q i PnR > f2(No) |, I*] < i, (1.43) 

where <, >pm ) is the real inner product on ^ 2 (No), P n is the projection on the set {n, n + 1, n + 2, . . . }, the 
operators A,B : £ 2 (Nq) — > ^ 2 (No) are defined by the kernels 

A(j,m) = (ip- 1 ) j+m+1 , B(m,k) = ip- m - k -i (1.44) 

and the functions Q, R € £ 2 (No) are given by 

Q(i) = (fV, m = ( — =/ bk+1 T^^ b )^r~ h - ( L45 ) 

\--z J\b\=i o-z 2mb 



On the other hand. 



n n (z) = z »e-X™Q°*»-»- k {l- < x _p nABPn P rJJ,P n V \z\ > 1, (1.46) 

where 

U(j) = (—^(-i) = I a-^^-I^, V{k) = (1.47) 
\z-- J j+1 J lal=1 z-aip{a)2ma 

Remark, (i) As 7T*(2;) is an entire function, the formula (|1.43|) also holds for a region of \z\ > 1 to which the 
right-hand-side of l|1.43fl is analytically continued, (ii) The conditions for <f> above can be weakened, but we 
do not discuss such an issue in this paper, (iii) Recall that Wicncr-Hopf factorizations of <j> are different by a 
factor of a multiplicative constant. However, since A and B have both factors -f^ -1 and tp, respectively, the 
operator AB is unaffected by a different choice of (j>+ and (f>- . The inner products in H1.43fl and i|1.46fl also 
remain the same even if ip is multiplied by a constant. Therefore, (| 1.43(1 and l(1.46JI does not depend on the 
choice of a Wicner-Hopf factorization of <f>. 

Remark. After this paper was completed, it turned out during a conversation with Andrei Martinez- 
Finkelstein that l|1.46fl appeared in [52] in a very different form. In the first formula of the equation (32) of 
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|22| . the authors found a series expansion for n(z). However one can check that the series is precisely what 
one would obtain once the Neuman series of the operator P n ABP n is taken in (|1.46(1 . The proof of [52] is 
based on a Riemann-Hilbert method. By turning the argument backward, it is also possible to prove the 
identity of Geromino-Case and Borodin-Okounkov using a Riemann-Hilbert method. This will appear in a 
future work. 

We regard 1)1. 43|) as an identity. We take a special choice of (f> and then take a limit of both sides of the 
identity i|1.43fl . A steepest-descent analysis shows that the right-hand-side converges to the formula i|1.15fl . 
On the other hand, a Riemann-Hilbert asymptotic analysis to the left-hand-side yields the Painleve formula 
(ll.2U|) . Hence the identity l|1.20fl follows from the identity between the orthogonal polynomials and their 
operator representation. 

This paper is organized as follows. In Section [21 we present the statistical and probabilistic models in 
which the distributions Fk(x;wi, . . . , Wk) appear. The proof of Theorem 1 1.1 1 is given in Section[3] Section^] 
proves Proposition II . 71 and Section proves Theorem 1 1.21 

Acknowledgments. The authors would like to thank P. Dcift, A. Its, A. Martincz-Finkelshtein, B. Simon 
and H. Widom for useful conversations and communications. This work was supported in part by NSF Grant 
#DMS-0350729 and the AMS Centennial Fellowship. 



2 Models 

We discuss several statistics and probability models in which F^'s appear. 



2.1 Non-null complex sample covariance matrices 

Let M > N > 1 be integers. Let y*x, ■ ■ ■ ,tjm be independent complex Gaussian N x 1 column vectors with 
mean /2 and population covariance E: the density of y\ is 

P[yi> (27r)^/2(det£)i/2 e [ > 

where * denotes the complex transpose. Denote by Y the sample mean Y := -^(j/i + • • • + Vm) an d by 
X = [y± — ?, . . . , i/M — Y] the (centered) M x TV sample matrix. Define the sample covariance matrix by 

S = y i xx *- ( 2 - 2 ) 

When the covariance matrix S is the identity matrix, the distribution of the eigenvalues of S is sometimes 
called the Laguerre unitary ensemble and is well-studied in the random matrix theory (see e.g. |12|). In 
particular, as M,N — * oo while M/N = r y 2 is in a compact subset of [l,oo), the largest eigenvalue A ma x 
satisfies the limit law (see e.g. [T3HT§] ') 

p((A max - (1 + 7- 1 ) 2 )) • (1+ 7 7)4 /3 M2/3 ^ x ) - *b(aO, (2.3) 
where Fg(x) is the Tracy- Widom distribution p. 3(1 . 
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Johnstone |2J| proposed the study of the so-called the 'spiked population model' where the covariance 
matrix E is a finite rank perturbation of the identity matrix. For possible applications of the spiked pop- 
ulation model in statistics, finance and telecommunications, see the references in and pp. For spiked 
population models, it is interesting to determine the effect of non-unit eigenvalues of the covariance matrix 
on the largest eigenvalue of the sample covariance matrix. For complex Gaussian samples, PQ determined the 
critical value of the non-unit covariance eigenvalue. When some of the non-unit eigenvalues of the covariance 
matrix are above the critical value, A max behaves differently from (|2.3|) . The function F k (x) is the limiting 
distribution of the A max when the largest eigenvalue of the covariance matrix is of multiplicity k and is equal 
to the critical value. 

Let £1 > • • ■ > £ r > be the non-unit eigenvalues of E where r is independent of M and N. 

Theorem 2.1 (Theorem 1.1 of pQ). As M,N — > oo such that M/N = j 2 lies in a compact subset of 
[l,oo), the following holds. 

(a) When 

£! = ■■■ = 4 = 1 + 7 _1 (2.4) 
for some < k < r, and t k +i, ■ ■ ■ ,&r are in a compact subset of (0, 1 + 7 _1 ), 

p((A max -(l + 7- 1 ) 2 ))- (1+ 7 7)4/3 M 2/3 <^ ^F k (x) (2.5) 

where F k (x) is defined in 

(b) When 

l\ = ■ ■ ■ = If. are in a compact subset of (1 + 7" 1 , oo) (2-6) 
for some 1 < k < r, and £k+i, ■ ■ ■ ,£ r ore in a compact subset of (0, £\), 




/ />2 7 -2 



P(^(A max - (h + ^-j)) ■ VM^jll - 1)2 < xj G k (x) (2.7) 

where G k (x) is the distribution of the largest eigenvalue of k x k Gaussian unitary ensemble. 

More detailed nature of the phase transition around the critical value 1 + 7 _1 was also studied in the 
same paper. 

Theorem 2.2 (Theorem 1.2 of [J). For some 1 < k < r, set 

i j = l + ^- {1+ £l Wj , j = l,2,...,k. (2.8) 



When Wi, . . . ,Wh are in a compact subset of R and tk+i, ■ ■ ■ ,£ r o,re in a compact subset of (0, 1 + 7 1 ), as 
M, N — > 00 while M /N = 7 2 is in a compact subset of [1, 00), 

p((A max -(l+7' 1 ) 2 )) ■ (1+ 7 7)4/3 ^ 2/3 <^ ^F k (x;w u ...,w k ) (2.9) 
where Fh(x; W\, . . . , Wk) is defined in 

It is transparent from this theorem that F k (x; w\, . . . , w k ) should be symmetric in wi, . . . , w k since re- 
labelling the eigenvalues docs not change the limit law. Further work on the eigenvalues of the spiked model 
can be found in |S] . 
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2.2 Last passage percolation and queues in tandem 

Suppose that to each lattice points € Z 2 , an independent random variable X(i,j) is associated. Let 

(1,1) / {N,M) denote the set of 'up/right paths' n = {(i k , jk)}^ 1 ' 1 where (i k+1 ,j k+1 ) - (i k ,jk) is 
cither (1,0) or (0,1), and = (1,1) and (in+m-i, .]n+m-i) = {N,M). Note that the cardinality of 

(1,1) /{N,M) is (™r 2 )- Set 

L(N,M):= max V X(i,j). (2.10) 

By interpreting X(i,j) as the (random) time spent to pass the site (i, j), L(N,M) is the Zasi passage time 
to travel from (1, 1) to (AT, M) along an admissible up/right path. 

Recall that the exponential random variable of mean m has the density function -^e~ x ^ m , x > 0. It is 
known that (see e.g. Proposition 6.1 of [P; we here scale X(i,j) of by M) when X(i,j) is an exponential 
random variable of mean £j (independent of j), nas the same distribution as the largest sample 

eigenvalue A max of complex Gaussian samples when the eigenvalues of the population covariance matrix 
£ are t\, . . . ,£n- Therefore for the last passage percolation model which have the identically distributed 
passage time for all but finitely many columns, Theorems 12.11 and Theorem 12 . 21 also hold with A max replaced 
by L ^} M ^ ■ In particular, Theorem 12.11 shows that as long as the site passage time on the distinguished 
columns have mean less than 1 + the last passage time has the same limit behavior as the case when 
all the sites are identically distributed. 

2.3 Queues in tandem 

Suppose that there are N servers and M customers. Initially all the customers are at the first server in a 
queue. Once a customer is served at a server, then (s)he moves to the queue of the next server and waits 
for his/her turn. The service time for the jth customer at the «th server is assumed to be a random variable 
X(i,j) and let D(N,M) be the departure time of all the customers from all the queues. It is well-known 
that D(N, M) has the same distribution as L(N, M) of the last passage percolation model (see e.g. |15|). 

In the queueing theory context, Theorem 12.11 determines the effect of a few slow servers to the total 
departure time. Suppose that X(i,j) is an exponential random variable of mean 1 for i = r + 1, N 
(independent of i) and of mean li for i = 1, . . . , r. In other words, the service times at the first r servers are 
distributed differently from those at the rest of the servers. When all of li are not so large, the departure time 
has the same limiting law as when all the service times are identically distributed, but the whole process 
slows down when some of the servers are sufficiently slow. Theorem 12.11 shows that the critical value is 
ti = 1 + 7 _1 . Note that due to a symmetry between servers and customers, the theorem also applies to slow 
customers. 



2.4 Totally asymmetric simple exclusion process 

The last passage percolation can also be interpreted as an interacting particle systems (see e.g. j20EI))- We 
will consider the totally asymmetric simple exclusion process. Let Xj(t) € Z, Xj(t), j = 1,2,..., t £ [0, 00), 
denote the location of the jth particle at time t. A particle can jump only to its right neighboring site after 
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random time if the site is not occupied. Let X(i,j) be independent random variables which represent the 
ith jumping time of the jth particle. We take the initial condition as Xj(0) = 1 — j, j = 1, 2, . . . . Then 
X(i 1 j) is the time it takes for the jth particle Xj to jump from the site i — j to i — j + 1. 

Let T(i, j) be the time it takes for the jth particle to arrive at the location i — j + Equivalently, T(i, j) 
is the time it takes for the jth particle to perform the first i jumps. Note that in order for the jth particle 
to jump from the site i — j to i — j + 1, the (j — l)th particle should be to the right of the site i — j + 1. 
Hence we find that 

T(i,j)=max{T(i-l,j),T(i,j-l)} + X(i,j), i,j>l, (2.11) 

where T(0, j) = T(i, 0) = 0, i, j > 1, by definition. A simple geometric consideration shows that last passage 
time L(i,j) satisfies exactly the same recurrence relation. Therefore T(i,j) is same as the last passage time 
in the sense of distribution. 

Let #(m,t) denote the number of particles to the right of the site m at time t. The flux F(m,t), the 
number of particles that have jumped cross the interval (m, m + 1) up to time t, is then F(m, t) = 4t={m, i) 
for m > 0, and F(m, t) = #(m, t) + m for m < 0. The event that #(m, t) > M is same as the event that the 
Mth particle is to the right of the site m at time t. This is again equal to the event that T(m + M, M) < t, 
and hence we find that P(#(m,t) > M) = P(T(m + M,M) < t) = P(L(m + M, M) < t). Therefore, 
Theorem 12.11 and Theorem 12.21 again apply to #([ut],t), and hence to F([ut],t). We state the results for 
#([ut],t) here. 

Traffic of slow start from stop 

Suppose that X(i,j) is an independent exponential random variable of mean li for i = 1, . . . ,r and of mean 
1 for i > r (independent of j). In other words, each particle jumps at rate j- for its first r jumps and then 
jumps at rate 1 afterwards. When t\ > ■ ■ ■ > l r , one can view it as a toy model for the following traffic 
situation: (infinite) cars in one lane, which were fully stopped at the red signal, speed up at the green signal 
until they finally reach the steady speed (after r 'jumps'). Set £ = max{^i, . . . ,£ r } and let k > 1 be the 
number of i^s equal to I. By re-interpreting Theorem 12. II a tedious but straightforward calculation shows 
the following results for — 1 < u < : 

( i o /i-it 2 \ 2/3 (f (-x), ue (1-7,01 n (-1,01 

HmP # {[ut],t )>t {l-u)H + x [^-) t 1/3 ) = { 1 ^ (2.12) 

V 4 V 4 J J | Ffc( _ x)j n = l - 1 e (-1,0], 

and for u £ (37, 1 - f) H (-1,0], 

f . . . . t-l-lu (£- lf/ 2 (l- 1 -lu) U9 \ , . 

& p (#a»*].*) >- -n^t + x w4 - 2 - w 2 ) = Gki ~ x) - 



(2.13) 



This shows that fast jumps do not affect the flux but slow jumps may change the flux. When r = (all cars 
jumping at the same rate), (|2.12|) was first obtained in ^5] for < u < 1. 
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Traffic with a few slow cars 

The exclusion process of the particles yields a dual process of the holes. As the particles jump to the right, 
the holes, the unoccupied sites, jump to the left. The leftmost hole jumps at rate j- since each particle jump 
at that rate at its first jump. Likewise, the second leftmost hole jumps at rate j- and so on. Hence this 
model can be thought of traffic model where there are a few cars of distinguished jump rates. Initially the 
holes are at the sites {1, 2,3,...}. As the number of holes Tt(m, t) on the left of the site to + 1 at time t 
satisfies Tt(m, t) = #(m, t) + m, (|2.12|) and (|2.13|1 imply the following results for — 1 < u < : 

/ 1 9 fl -u 2 \ 2/3 w ,\ {f (-x), u € (1 - |,0] n (-1,01 

lim F[H([ut],t) > -(l + u) 2 t + x( — - — ) t 1/3 = S (2.14) 



'-* 00 V Iv.r ~\ u = l-| G (-1,0] 




and for u e 1 - |) n (-1,0] 



,//r i \ £-l-£u + e 2 u (£- lf/ 2 U- l-iu) 1/2 \ . ,„„„, 

«([«*],*)> ^ * + * V(^-2-M 1/2 J fc( ~ X) ' ( } 



i — >oo 



The full case of —1 < u < 1 and also correlation functions of various locations for both of the above 
traffic models will be discussed in a forthcoming paper. 



3 Proof of Theorem 11.11 

We prove Theorem II .11 in this section. 

Since both sides of (|1.16fl are analytic in each Wj, the case when some of Wj's coincide follows from 
analytic continuation of the case when all Wj's are distinct. Hence we assume in this section that all Wj's 
are distinct. We need to prove that 



<(^> mn - < 1 . . .,W m ),t^(Wl, ■ ■ ■ , W n -l) >(x,oo) 

dct {{w m + D x ) n - 1 f{x,w m ) 



det| S mn - < n 

1 < m , n < k 

1 , ( , ^ 1 „ x\ (3- 1 ) 

A^-i- I (... i t-» \n.— \ _(*/___ . \ t v 7 



(io„ - w m ) 



Km,n<k 



l<m<n<fc 

Notational Remark. In the below, we sometimes have a product of empty indices. For instance when 

71-1 

n = 1, the product (w a — «;,) in 13.11(1 has no indices. In such cases, we interpret the product as 1. 

CJ=1 

Let A x : L 2 ((0, oo)) — > L 2 ((Q, co)) be the operator with kernel 

A^u) = A(m + x,w + x). (3.2) 

Set 

S m (u) = s (m) (u + x) = s {m) (u + x;w u ...,w m ) (3.3) 

and set 

T m (v) = & n \v +x) = & n \v +x;w 1 ,.. .,w m -i). (3.4) 



13 



Then 

(t^a:»"" , )<««>-(t J a: s - 

and the matrix on the left-hand-side of i|3.1[) is 

( s v-<T^ s i> T i>)i<i, j <k 

where <, >=<, >(o,oo) is the rea l inner product in £ 2 ((0, oo)). 
Since 

rn m 1 m 1 

I r . = y - | r - , 

- 1 - 1 - Wj + ia ^— ' Wj + ia uup — Wj' 

j=i j=i J i=\ ■> 



we find 



where 



m r m 



' 3 = 1 '-£=1 



E w (u) = E w {u;x) 



1 — A, 



Cu, (u), C7«,(u) := C w (u + x) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(recall (I1.14|) for the definition of C w ). For later use, we note that f(x,w) defined in i|1.15fl satisfies that 

f(x,w) = 1- < E w ,Ti > . (3.10) 
Now we invert the relation Ij3.8|l . For 1 < m < k, l|3.8|l is a system of k linear equations for E Wj , 1 < j < k. 
Lemma 3.1. The equation ()3.8|) for E Wj has the solution given by 

1 



3 rn-1 



E Wj (u) =^ Y[(Wa ~ Wj) 



n—1 L a— 1 



l - . 



s„ (it). 



Proof. Consider the function 



m 

F(z) : = 1] — 



— z 



(3.11) 



(3.12) 



Integrating over a circle of radius R, and then taking R — > oo, we find that the sum of residues of F is equal 
to when m — n > 1 and is equal to 1 when m = n. On the other hand, by directly computation, the residue 



of F at z = Wj is 



Hence we obtain the identity 



f X Wi - Wj 

£—n J 



mm ^ 

Z] TT = Smn, 

* — ' ■ L1 - Wl — Wj 

j=n l=n J 



m > n. 



(3.13) 



(3.14) 



Now as all w^s arc distinct, the determinant of the matrix for the linear equation l|3.8|) is Y[i<e< m <k( w £ — 
Wm) -1 , which is non-zero. Hence there is a unique solution E Wj for l|3.8|) . We should check that (|3.11|) solves 
(|3.8|) . But this follows by inserting (|3.11|) into the right-hand-side of 13. 8f) . changing the order of summations, 
and then using 1|3.14[) . □ 
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From l |3.11J k we obtain for each 1 < i, j < k, 

i 7i — l 



-S n ,Tj >) — Fij— < E Wi , Tj > 



n—1 a—1 



where 



j'-i 

a=l 



(3.15) 



(3.16) 



Note that — when i < j . Now we perform row operations of the matrix (<%— < j-r^-Si, Tj >) 1<i j <k 
using (|3.15|l that replaces the zth row by a linear combination of the first i rows to find that 



k i-l 



L 3 > )kxk 



det(F J -<E m ,T J >) kxk . 



= 1 a=l 



(3.17) 



n (^ _w j) 

l<i<j<fc 



Note that the when j = 1, = 1, (see the Notational Remark above) and hence the first column of the 
matrix (Fij— < E Wi ,Tj >) kxk consists of the functions (see (13.101) ) 



1- < E Wi ,Ti >= f(x,Wi). 
For example, when k = 3, the determinant on the right-hand-side of (|3.17l) is 



(3.18) 



det 



f(x,wi) - < E Wl ,T 2 > 

f(x,w 2 ) (wi - w 2 )- < E W2 ,T 2 > 



— < E W1 , T s > 

— < E,., n , T 3 > 



\ 



(3.19) 



\f{x,w 3 ) (wi - w 3 )- < E W3 ,T 2 > (to! -w 3 )(w 2 - w 3 )- < E W3 ,T 3 >J 
From the definition l|3.4|l of Tj and the definition (|1.10|> of tj , we have 

Tj = Wj-\Tj_i - D x Tj-i, j > 2. 

Set M(°) be the matrix 

m(°) : (/■;, f„..t, • ),..,.,. 

Let be the matrix defined by 



(3.20) 
(3.21) 



(\ -w x 
1 -w 2 

1 -w 3 



M (D := M (o) 



\ 



1 -w k -i 
1 / 



(3.22) 
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whose determinant is same as the determinant of M^°\ Using the relation l|3.20fl . the entries of 



V )\<i,j<k 



are given by 



M 



(i) _ J f( x , w i), 



where 



F^ — < E Wi ,D x Tj-i >), 



F iP : " w 3-l F i,j-l ~ F ij- 



3 = 1 

2 < j < k, 



Now define a new matrix M^> — (My ) 1<;i <k by 



y )\<i,i<k 
(l 

1 -K/i 

1 -w 2 



M< 2 > := M« 



V 



1 -Wk-2 

o i y 



Using the relation 

that follows from i|3.2U[l for 3 < j < k, we find that 

f{x,Wi), j = l 

V = \-{F£ ) -<E m ,D x T 1 >), j=2 

,(2) 



F£> - < E Wi , D 2 x Tj- 2 >, 3 < i < fc, 



where 



(2) 



(i) _ F (X) 

3-1 *J ' 



Wj- 2 FM_, - F 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



1G 



Continuing in a similar way, we eventually define 



/I 





M (fe-D := M (o) 



-Wl 

1 







1 





-w 3 



(I 
1 




/l 
10 
1 



-Wk-l 

1 

\ 



1 





-U>2 



V 



1 




1 / 



Using the fact that for all t > 1, j > 2, 
we find that 



71/ 



= (-ir 1 ^?"^- < E^D^T, >), 1 < i, j < fc, 



where for each 1 < i < k, is inductively defined by the relation 



F, 



(0 



and (see (|3.16[) ') 



w . eF v- 1 ) _ F v-v 



1 < ^ < j < fc 



4 := F« = K - 1 < j < fc. 



Recall that i 7 ^ = 1, and hence ^ = 1— < E Wi ,T\ >= f(x,Wi). 

Lemma 3.2. The solution F^ to the recurrence relation <|3.32[) and <|3.33[) is 

3 -l-l 



F^=wl 



[ (w a - Wi ), i<e<j<k. 



-Wk-2 
1 / 



(3.29) 

(3.30) 
(3.31) 

(3.32) 

(3.33) 



(3.34) 



Proof. This follows easily from an induction in I. Here, as mentioned in the Notational Remark above, when 
j = I + 1, we understand that the product n^iC^a — w i) = 1- '-' 
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Therefore F^~ 1] = w 3 ' 1 , and as det(M^) = de^M^" 1 )), we find from (|3~T7|) . j3~2TJ) and (|3~3T|) that 



< T^A^ >W - jj ^ dCt(M) ' (3.35) 

l<i<j<fe 

where [fc/2] denotes the largest integer smaller than or equal to fc/2, and the k x fc matrix M = {Mij) 1<i - <k 
is given by 

= wf 1 - < E w ^D j x ' 1 T 1 > . (3.36) 
As ni< m <„< fc (-l) = (-l) lk/2] , this is equal to 



det(%- < ^—^-^.Tj- >) 1 < lJ < fe 



det(M) 



n k - ( 3 - 3? ) 

l<i<j<fc 



Now we will show that det(M) is equal to the determinant on the right-hand-side of (|3.1J) . For this 
purpose, we use the following result. 

Lemma 3.3. For I > 0, there are smooth functions F^ a (x), a = 0, 1, 1, smc/i t/iat 

(w + IkjVfou;)^- < ^,^Tx >-£i^ o (a0(u; + i? x ) a /(z,«0- (3-38) 

a=0 

Proof. From the definition of E w , 

d x £„, = d x (—L-cS) = Y-^( D * k ^Y~hk^ w + T^T Dx ^ w ' (3 ' 39) 

It is direct to check that (D x A x )(u,v) = — Ai(x + u) Ai(x + v). Hence D X A X = —T\ ® T\. On the other 
hand, from the definition i|1.14[) of C w , D X C W =T\ — wC w . Hence we find 

D X E W = -—L—Ti ® T x - — C w + Y^r-iTi - wC w ) 

= _ < Tu rhr x dw > t 1 ^/ 1 + i^r x Tl w T^r x dw ' 

which implies that 

(w + D W )E W = f(x, w )y^~~A~ Ti (3 ' 41) 
Now we use an induction in I to prove (|3.38(l . When I = 0, by definition (|1.15() of /, i|3.38|l holds. Now 
suppose that l|3.38|l holds true for some i > 0. Then using the general identities (w + D x ) < h,g >=< 
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(w + D x )h,g > + < h,D x g > and (w + D x )(hg) = (D x h)g + h({w + D x )g) , 

(w + D x ) e+1 f(x,w) 

= (w + D x )[(w + D x ) e f{x lW )] 

1-1 

= {w + D x ) [w l - < E w ,D e x T 1 >-J2 Ft, a {x){w + D x ) a f(x, w)] 

a=0 

= w e+1 - < (w + D X )E W , D e x Ti >-< E w ,Di +1 T! > 

- W (D x F i!a (x)) (w + D x ) a f(x, w) + F t , a (x){w + D x ) a+1 f(x, w) 



(3.42) 



= w e+1 - f{x, w) < j^T x ,DiT x >-< E w ,D i + 1 T 1 > 
e-i f 

- J2\ (D x F i>a (x)) {w + D x ) a f(x, w) + F^ a (x)(w + D x ) a+1 f{x, w) 

a=0 ^ 

where (|3.41|l is applied in the last step. Therefore we find that l|3.38[) holds true for £ + 1 with the functions 

D w Ft,o(x)+ < t^Ti, D l x T x >, a = 0, 

{ D x F La {x) + F tta -x, l<a<£-l, (3.43) 

Ft,e-i, a = £, 
where Fq^ := 0. □ 

By applying l|3.38[) repeatedly to (w + D x ) a f (x; w) inside the summation on the right-hand-side of 1|3.38[) . 
for £ > 0, there are smooth functions Gi a , a = 0, 1, ...,£— 1 such that 



w £ - < E w ,D x Tt >= (w + D x )'f{x- w)+J2 G e ,a{x)(™ a - < E w ,D a x T x >). (3.44) 

a=0 

Therefore for any 1 < i < fc, 

Mij = wj- 1 - < E^Di- 1 ^ >= (wi + D x ) j ~ 1 f(x,w i ) + ^2 Gj-i,a-i(x)M ia . (3.45) 

a=l 

In other words, the jth column vector in the matrix M is equal to a linear combination of the first, second, 

j — 1th column vectors plus the vector ((wi + D x y^ 1 f(x,wi), • • • , (wk + D x Y~ 1 f{x,Wkj) T ■ Hence by 
applying proper column operations, we find 

det(M) = det(K +D x y- 1 f(x,w i )) 1 <. t (3.46) 

This, together with i|3.37[l . implies that the left-hand-side of i|3.1|) is equal to 



1 



det( K (w m +D x ) n - 1 f(x,w m )) i 



(w n -w m ) ^ ' m ^i<m,n<k (3.47) 



l<m<n<fc 

This is the right-hand-side of (|3.1() and Theorem 1 1.1 1 is proved. 
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4 Proof Proposition 11.71 



Let T n (<p) = (ipi-j)o<i,j< n -i be the Toeplitz matrix of the symbol ip where tfk denotes the Fourier coefficients 
of ip. Let D n (ip) = detT n (p) be the Toeplitz determinant. We will use the following identity ^][7] between 
a Toeplitz determinant and the Frcdholm determinant of an operator on a discrete set: for all n > 1, 

Gwm =det{1 - PnABPn) (41) 

with 

( 00 ^ 

(4.2) 



X>(log^) fc (log^)- 
fc=i 



where the operators A, B are defined by the kernels (|f .44f) . This identity holds, for example, for complex- 
valued analytic functions tp with zero winding number which has a Weiner-Hopf factorization. See [§] for 
the minimal condition on ip for which the identity holds. 

It is well-known that 7r* has the multi- integral expression (see e.g. |26| ) 

<{z)=^- j .. [ n (i-^na-^n^)^. m 

Vn(<P) J\ Z1 \ = 1 J\ Zn \=! x < j<k < n j=1 j=1 

By using the multi-integral formula of a Toeplitz determinant, 1|4.3[) can be written as 



where the new symbol (f> z is 



z 



See |2] for a use of the identity l|4.4|l in random matrix theory. Using (|4.1fl for D n (<f> z ) and D n ((/>), we find 
that 

G(dp z ) n E{<j) z ) det(l — P n A z B z P n ) 



G{4>) n E{(p) det(l - P n ABP n ) 



(4.6) 



G((f>) n E(4>) V I -PnABP n 

where the operators A z , B z : £ 2 (Nq) — > ^ 2 (No) arc defined by the kernels 

A z (j, m) = (ip~ 1 ) j+m+1 , B z (m,k) = (ip z )- m - k -i. (4.7) 

where ip g = (0 Z ) + /(0 Z )_. 

When \z\ < 1, from (|4.2|l . it is easy to check that 

G{<j>)»E(4>) 

Now we consider A Z B Z — AB. As Z has the Wicner-Hopf factorization Z = (0 Z ) + (0 Z )_ where (^> z )+(w) 
4>+{w) and (0 z )_(u>) = (l — ^)</>_(w), we find ip z (w) = ■jj^j^j'4>{w) . Therefore, 



(4.8) 



(A z B z )(j, k) - (AB)(j,k) = ^2(lp z 1 )j+,n+l(lp Z )-rn-k-l - {ip 1 )j+rn+l^-rn-k-l 

= f f a -j- m -x b m+k+i( lzjk da db 

±^J\ a \=iJ\b\=i U-f / i>(a) 2nia 2nib 
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Note that 4>+{z) = 4>(z)/(f>-(z) is analytic in a region outside the unit circle as <fr(z) is assumed to be analytic 
in a neighborhood of the unit circle and <f>-(z) is analytic outside the unit circle. Hence <j>+(z) is analytic in 
a neighborhood of the unit circle. Similarly, (/>_ (z) is analytic in a neighborhood of the unit circle and so is 
ip(z). Therefore the contours \a\ = 1 and |6| = 1 can be deformed so that \a\ > \b\. Hence, 



(A z B z )(j, k) - (AB)(j, k)= a-i-'b^ 1 

l\a\>\b\ 



r oc / 7 




E 

m=0 



|o|>|6| b - z Via) {2m) 2 ab 



z{a — b) ip(b) dadb 

a(b — z) ib(a) (2'Ki) 2 ab 

' W V ' (4.10) 



where Q and R arc defined by l|1.45f) . This implies that A Z B Z is a rank 1 perturbation of AB and we find 

<(z) = e -ES=iO»K*)*« fc .det(l - 1 _p nABPn P nQ®RP^j 



(4.11) 



which completes the proof of (|1.43|) . 

Proof for l|1.46|l is similar by noting that 

*n(z) = ^T^, <f(w):=(z-w)<Kw). (4.12) 

5 Proof of Theorem 11.21 

We apply Proposition II .71 to the function 

:= e*<*+* -1 > (5.1) 

for positive number t. Then l|1.43|l becomes 

e tz <(z) = 1- < —A x —P n Q, Pn R >, 2(Nq) (5.2) 

with 

^(s) = e*( z -^). (5.3) 

It is easy to check that the inner product on the right-hand-side is unchanged when the functions A(j, m), 
B(m, k), Q(j) and R(k) are replaced by 

(-iy +m A(j,m), (-l) m+k B(m,k), (-l) j Q(j), (-l) k R(k), (5.4) 

respectively. We will denote these new functions by the same notations A,B,Q,R. 
We will take the limit t — ► oo in both sides of the identity 15.2fl with the scaling 

n = [2t + xt 1/3 ], z = -\ + ™ (5.5) 

for a fixed real number x and a complex number w, where [a] denotes the largest integer smaller than or equal 
to a. We will see that under this scaling limit, the right-hand-side of the identity (|5.2|l becomes f(x, w) given 
in H1.15[) and the left-hand-side becomes l|1.20fl . thereby yielding the desired Painleve formula of f(x,w). 
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Indeed, the limit of e tz 7r* (z) is obtained in |J and also [3]. The paper 0] studied the asymptotic behavior 
of the longest increasing subsequences of certain symmetrized versions of permutations and the asymptotic 
analysis of e* z 7r*(z) was a technical part of the paper. The paper [3] on the other hand studied a last 
passage percolation model, which is different from the one discussed in Section|21 From (5.22) and (5.26) of 
Proposition 5.4 and Corollary 5.5 of 0], we find that 

I M 2 2(-Uw,x), w>0 
y-M 2 i{-\iw;x)e^ w ~ xw , w < 0. 

This result actually motivated us to use the function Q5.1[l . 
On the other hand, it is known that [SI [20] ( see a ^ so |29| ~) 

P n ABP n -> A, (5.7) 

in trace norm for any fixed real number x where A x is the Airy operator defined in 1(1. 2|) . This limit was 
studied in the papers [HI 123 GUI in the context of the longest increasing subsequences and the Plancherel 
measure on partitions. Therefore, the only remaining part is the asymptotic analysis of Q and R. These 
can be done by a standard steepest-descent analysis. Similar analysis appeared in several places (see e.g. 
[HlEin) an d we only sketch basic ideas. 

We only consider R since the analysis of Q is similar. Note that the integral formula l|1.45f) . which was 
originally defined for \z\ < 1, can be analytically continued for all complex numbers z by deforming the 
contour so that z lies inside the contour. To compute the limit of the right-hand-side of (|5.2(l . we need the 
limit of R([2t + yt 1 ^ 3 ]) with certain uniformity for y G [x, oo) to ensure the convergence of the inner product. 
As one can check from the analysis, it is reasonable to think that R([2t + yt 1 / 3 ]) is close to R(2t + yt 1 ^ 3 ) and 
we will compute the limit of the later. See \n^, f° r example, for a discussion of this type. Now from l|5.3|) 
and 15.5(1 

«< 2t + *"'-> - h I t H-i^./3) '- t » 2W " % " t ' t " , ' fa < 5 - 8 > 

where the contour is modified to go from oo + iO to oo — iO enclosing the origin and the point — 1 + wt^ 1 / 3 . 
Here the function (— b) 2t+yt ' denotes the principal branch. Note that the integrand is of the form 

- 1 + wt ~ 1/3 . { - br ^ e t m (59) 



b - (-1 + wt-V 3 ) 
where 

/(6) = 21og(-6) + &-fe- 1 . (5.10) 

The function f(b) has the double critical point at b = -1, and /(-l) = f'(-l) = f"(-l) = and / (3) (-l) = 
2. Approximately the steepest-descent contour passing the critical point b = —1 is, in a neighborhood of 
size, say e > 0, of b = — 1, the union of the line from — 1 + ee 7 ™/ 3 to —1 and its complex conjugate. As the 
pole — 1 + w 1 / 3 lies to the right of b = — 1 when w > 0, one can also check that it is possible to deform the 
original contour to the steepest-descent contour when w > 0. When w < 0, we can modify the contour to 
the union of { — 1 + xe 7 ' 1 / 3 : 2\w\t~ 1 / 3 < x < e}, {2\w\t~ 1 / 3 e ie : n/3 < 8 < n} and their complex conjugate so 
that the pole — 1 + wt^ 1 / 3 still lies on the right of the contour but the contour 'essentially' passes the point 
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b = — 1. See PP where a similar modification of the contour was used for a steepest-descent analysis. Note 
that in the neighborhood of b = — 1, the contour is oriented from — 1 + ee 7 ™/ 3 to — 1 + ee~ 7 ™ / ' 3 . 

From the standard steepest-descent method, the integral is asymptotic to the integral over the part of 
the contour in the e neighborhood of a = — 1. The approximation f(b) ~ gr/ ( — + I) 3 = — §(& + I) 3 
suggests the change of variables zt 1 / 3 (6 + 1) = s, which implies that 



~ — / — - — e i!/s+i 3 s3 ds 
2?r 



(5.11) 



where the contour is from ooe 57 ™'' 6 to ooe 7 ™'' 6 such that the pole s — iw is above the contour. Hence wc find 
that R(2t + yt 1 ^ 3 ) ~ —C w (y). Similar calculation shows that Q(2t + yt 1 ^ 3 ) ~ — Ai(y). This argument can 
be made rigorous with uniform error bound for y (see e.g. pQ for a similar calculation). Therefore, by noting 
that 1 — A x is a self-adjoint operator, one finds that the right-hand-side of l|5.2|l converges to (|1.15|) . Thus 
we obtain the identity 

1 I M22(—hiw;x), w>0 

l-<^^-C w ,Ai> LH(x>oo)) ={ 2 i 3 (5.12) 

1 - A * [-M2i(-±iw]x)e3 w '- xw , w<0. 

The proof of Theorem 1 1.21 is complete. 
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